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FORMULAS  FOR  THE  PRESSURE  AND 
BULK  MODULUS  IN  UNIAXIAL  STRAIN 

Mike  Scheidler 

U$  Army  Research  Laboratory,  APG,  MD  21005-5066 

For  an  isotropic  elastic  solid,  the  pressure  p  =  pu(p)  in  a  state  of  uniaxial  strain  at  density  p 
generally  differs  from  the  pressure  p  =  ph(p)  in  a  state  of  hydrostatic  stress  at  the  same  density. 
Several  researchers  have  used  pressure/shear  (or  oblique  plate  impact)  tests  to  determine  pu  and 
the  corresponding  uniaxial  bulk  modulus  =  pdp^/dp,  The  pressure/shear  tests  yield  uniaxial 
longitudinal  and  shear  moduli,  Lu  and  GUi  as  a  function  p,  A  common  procedure  is  to  integrate  the 
approximate  relation  Kv  &  -  | Gu  to  obtain  the  pressure-density  relation  p  —  pu(p)  in  uniaxial 

strain.  It  is  shown  here  that  the  integration  of  this  approximate  relation  between  the  moduli  can 
be  avoided  altogether  by  utilizing  the  exact  formula  pu  =  -  | ((p/pa)2  —  l)Gu,  where  01  denotes 

the  longitudinal  stress  (pos.  in  compression).  The  bulk  modulus  is  computed  exactly  from  this 
formula,  and  the  error  in  approximating  it  by  Lu  —  is  determined. 


1.  INTRODUCTION 
We  consider  only  isotropic  elastic  materials, 
and  for  simplicity  thermal  effects  are  neglected 
until  §5.  Under  these  conditions,  the  pressure  is 
typically  assumed  to  be  a  function  of  density  only. 
However,  nonlinear  elasticity  theory  predicts  that 
the  pressure  also  depends  on  the  shear  strain,  al¬ 
though  isotropy  implies  this  effect  is  necessarily 
of  second  order;  cf.  Scheidler  (1),  In  §4  we  derive 
exact  formulas  for  the  pressure  and  bulk  modu¬ 
lus  in  a  state  of  uniaxial  strain.  The  effect  of 
shear  strain  can  be  seen  by  comparing  these  re¬ 
sults  with  the  corresponding  relations  for  a  state 
of  hydrostatic  stress  (§3),  Our  results  are  based 
on  exact  formulas  for  the  speeds  of  acceleration 
waves  (§2).  Applications  to  the  analysis  of  data 
from  pressure /shear  tests  are  discussed  in  §5. 

2.  ACCELERATION  WAVE  SPEEDS 
Let  F  denote  the  deformation  gradient  relative 
to  the  undeformed  and  unstressed  state.  The  left 
Cauchy- Green  tensor  B  =  FFT  has  principal  val¬ 
ues  hi  =  V,  where  As  are  the  principal  stretches, 
and 

detF  =  —  A1A2A3  =  — ,  p  =  — ,  (1) 

p  Po 


where  p  and  p 0  denote  the  densities  in  the  de¬ 
formed  and  undeformed  state.  The  principal  axes 
of  B  are  the  principal  axes  of  strain  in  the  de¬ 
formed  state.  Since  the  material  is  isotropic  and 
elastic,  these  axes  are  also  the  principal  axes  of 
the  Cauchy  stress  tensor  T,  and  T  is  an  isotropic 
function  of  B.  This  implies  that  there  is  a  single 
function  t  such  that  the  principal  stresses  i*  are 
given  by 

=  (?) 


for  any  permutation  i,jt  h  of  1,2,3;  cf.  Truesdell 
&  Noll  (2).  It  follows  that  the  pressure 

p  =  —  |trT  =  —  |(*i  +  h  +  *s)  (3) 


is  a  symmetric  function  of  6it  63163-  Analogous 
results  hold  in  terms  of  the  principal  stretches  Ai 
or  in  terms  of  various  principal  strain  measures, 
e.g.,  A,  -  1.  tfn  -  1).  A(1  -  1  /Ai),  or  In  A,  . 

The  speed  Ui  of  a  longitudinal  acceleration 
wave  propagating  along  the  ith  principal  axis  of 
strain  is  given  by 


pU?  =  lbi 


dU_ 

dbi 


.dtj  dti 

*dX d  In  A, 


W 


The  speed  V\j  of  a  transverse  or  shear  accelera¬ 
tion  wave  propagating  along  the  ith  principal  axis 
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of  strain  with  jump  in  acceleration  parallel  to  the 
jth  principal  axis  (j  ^  i)  is  given  by 


=  if  M&j-  (5) 

All  quantities  in  (4)  and  (5)  are  evaluated  at  the 
wave  front.  These  wave  speeds  are  in  the  de¬ 
formed  material  (i.e.,  Eulerian);  the  correspond¬ 
ing  Lag  ran  gi  an  wave  speeds  arc  obtained  by  di¬ 
viding  by  Xi.  Proofs  of  (4)  and  (5)  can  be  found 
in  Truesdell  k  Noll  (2)  and  Wang  k  Truesdell 
(3).  These  formulae  do  not  require  that  the  region 
ahead  of  the  wave  be  at  rest  or  in  a  homogeneous 
state  of  strain.  However,  when  these  conditions 
are  satisfied,  (4)  and  (5)  also  apply  to  the  speeds 
of  plane  infinitesimal  sinusoidal  waves;  cf,  Trues¬ 
dell  k  Noll  (2). 


3,  HYDROSTATIC  STRESS 

For  a  purely  dilatational  deformation,  we  have 

i>i  =  p~2fz  and  £*  =  —p  (i  —  1,2,3),  (6) 

In  this  hydrostatic  stress  state  every  axis  is  a  prin¬ 
cipal  axis  of  stress  and  strain,  and  the  pressure 
p  is  a  function  ph  of  p  or  p.  Here  and  below,  an 
“h”  subscript  denotes  the  hydrostatic  stress  state. 
From  (2),  (4),  (5)i ,  and  (6),  it  follows  that  for  a 
given  density  p  there  is  a  single  longitudinal  wave 
speed  Ui  —  Ui^h  and  a  single  shear  wave  speed 

Uij  =  Us, ht  an(J  that 


dph 

dp 


uL>h  ~  * 


cf.  Wang  k  Truesdell  (3).  A  different  proof  of  this 
well-known  result  is  given  by  Truesdell  k  Noll  (2). 
We  assume  ph  is  a  strictly  increasing  function  of 
p.  Then  (7)  implies  the  longitudinal  wave  speed 
is  greater  than  the  shear  wave  speed.  With  the 
longitudinal,  shear,  and  bulk  moduli  defined  by 


(7)  implies  the  well-known  relation 

Kh=Lh-  f  Gh  .  (10) 


We  use  a  zero  subscript  to  denote  functions  eval¬ 
uated  at  the  undeformed  state  where  A,  —  bi  — 
p  —  1;  in  particular,  Kq  —  La  —  By  (9), 

M/ffoMp-n  +  ietf-lJ*.  (11) 

where  the  dimensionless  constant  oo  is  given  by 


do  as 


1  <l2ph 


2A'0  dp 1 


4(tH  « 


4.  UNIAXIAL  STRAIN 

For  a  state  of  uniaxial  strain  along  the  l-axis, 

Ai  —  y/b\  —  1  /pi  X2  =  X$  —  Aj  —  £3  =  1  (13) 

and  (2)  implies  —  £3.  The  principal  stresses 
^  are  positive  in  tension;  if  =  —  U  then  cq 
is  positive  in  compression.  We  use  a  Hti"  sub¬ 
script  to  denote  uniaxial  strain  and  consider  only 
waves  propagating  along  the  l-axis  into  a  uniax- 
ially  strained  material.  The  Eulerian  wave  speed 
Ul, u  —  Ui  of  a  longitudinal  acceleration  wave  is 
given  by  (4)  with  i  ==  1,  and  hy  (13)  we  also  have 

LuSpV^=pd^  =  P^..  (14) 

It  follows  that  a  longitudinal  acceleration  wave 
can  propagate  only  if  dtrifdp  >  G,  Le.,  if  is 
a  strictly  increasing  function  of  p,  which  we  now 
assume.  By  (13),  the  material  is  strained  iffp  ^  1 
iff  61  ^  63.  In  this  case  (5)s  and  (13)  imply  the 
following  formulas  for  the  Eulerian  speed  U$tii  = 
f/u  of  a  transverse  or  shear  acceleration  wave: 


Gh  =  p  USu  —  hx 


£1  —  h  *i  —  £2 


i  — 

“  p-1  ~ 
where  r  is  the  shear  stress: 


 I -p 

2  r 

P2-V 


(15) 


L>,=pV&  Gh  =  PUs\  (8) 


r  =  i(<ri  —  cr2)  =  (16) 


The  Lagrangian  wave  speeds  are  p  Ui^u  and  pUs^u- 
If  p  >  1  the  material  is  in  compression,  and  (15) 
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implies  that  a  shear  acceleration  wave  can  prop¬ 
agate  only  if  <Ti  >  tr2.  When  p  —  1,  the  re¬ 
sults  of  the  previous  section  apply,  and  we  have 
Gu jo  "  Go  and  Lu\q  — 

From  (15}  we  have  the  following  fundamental 
formula  for  the  shear  stress  in  uniaxial  strain: 

T=±(p2-l)Gu.  (IT) 

Since  t2  —  *3,  (3)  and  (16)  imply  the  follow¬ 
ing  well-known  relation  between  the  (compres¬ 
sive)  longitudinal  stress  <7^  the  shear  stress  r, 
and  the  pressure  p*.  in  uniaxial  strain: 

P«  =  <T\  -  f  T  .  (18) 

On  substituting  (17)  into  (18),  we  obtain  the  fol¬ 
lowing  fundamental  formula  for  pu: 


pu=ff!-|(p2 -!)(?-.  (19) 

We  define  the  bulk  modulus  Ku  in  uniaxial  strain 
by 

v-  —  dPu  _  -  _  j  ^Pv  /9n\ 

h'  =  p7i  =  plp-L'‘d71  (20) 

Then  from  (19)  and  (14),  we  obtain 


The  results  up  to  this  point  are  exact.  We  now 
consider  some  useful  approximate  relations.  FYom 
(20)  we  have 

Pu/A'o  «  (p  —  1)  +  bo(p  —  l)3,  (24) 


where  the  dimensionless  constant  b0  is  given  by 


_L£Y“| 

2Ao  dp*  |0  “  2  V  dpu  |„  . 


(25) 


For  use  in  (24)“{25),  note  that  (21)  implies 

d A  u  Lq  (  dLx,  I  8  dGu  |  \  8  Go 

dpv  K0\  dffX  |a  3  d<n  3  AV 

(26) 

From  (11 )“( 12)  and  (24)— (25) ,  we  see  that 

Pu  »  Ph  +  KaCo(p  -  l)3  .  (27) 

^  P^Pft  ^  £o(  *_  _  (28) 

Pu  ph 

where  the  dimensionless  constant  Co  is  given  by 


Co  —  fro  —  do  = 


1  (  dKu  I 

2  \  dpu  |Q 


(29) 


A«  =  £«-|p3G.-|(P3-1)P^ 

=  (L0-fG„)-|(p2-l)fi„,  (21) 

Hu  =  2 Gu  +  p™  =  ~(p3G„)  ■  (22) 

dp  p  dp 

At  p  —  1,  (21)  reduces  to  Ku\q  =  Lq—  |Gq  =  A'o- 
From  (21)  and  (22)  it  follows  that  A'„  =  Lu  -  |GU 
for  all  p  iff  Hu  =  0  iff  <7^  =  Go/p2,  but  there  is  no 
reason  to  expect  such  dependence  in  general,  and 
thus  no  reason  to  expect  that  Ar„  =  Lv  -  |Gt, 
except  in  the  limit  of  aero  strain.  Of  course, 
by  analogy  with  (10)  we  could  have  defined  A„ 
to  be  Lu  —  |Gti,  but  then  (20)  would  not  hold. 
From  (21 )  we  see  that  for  a  state  of  compression, 
Ku  <  Lu  —  iGu  if  Hu  >  0,  and  Ku  >  Lu-  %GU  if 
<  0.  We  assume  that  pu  is  a  strictly  increas¬ 
ing  function  of  p.  Then  any  function  of  p  may 
also  be  regarded  as  a  function  of  &i  or  pUi  and  by 
(14)  and  (20)  we  have 


_  cf  _  d  _  r  d 
dp  **  dp  d<T\  dpu 


(23) 


On  comparing  (28)  with  equation  (4.6)  in  Schei- 
dler  (1),  we  find  that  cq  is  also  given  by 


_  2  fdGh  1  Go  \ 
"  3  IdpJ*  An/ 


(30) 


6.  DISCUSSION 

The  longitudinal  stress  ci  as  a  function  of  p 
in  uniaxial  strain  can  be  obtained  from  normal 
plate  impact  tests.  Then  the  relation  (18)  (which 
does  not  rely  on  the  assumption  that  the  response 
is  elastic)  is  typically  used  to  determine  the  pres¬ 
sure  pu  in  uniaxial  strain  given  some  assumptions 
on  the  shear  stress  rt  or  to  determine  r  given 
some  assumptions  on  pu^  It  is  often  assumed  that 
pu(p)  is  equal  to  the  pressure  ph(p)  in  a  state  of 
hydrostatic  stress  at  density  p  (or  to  some  ap¬ 
propriate  modification  of  p^  to  include  thermal 
effects  in  the  shocked  state).  Such  an  approxima¬ 
tion  neglects  the  effects  of  shear  strain  (or  shear 
stress)  on  pu .  That  this  effect  may  be  significant 
in  ceramics,  geologic  materials,  and  polymers  has 
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been  emphasized  by  Gupta  (4)  and  Conner  (5). 
These  materials  can  sustain  relatively  large  elastic 
shear  strains  (compared  to  metals),  although  for 
polymers  viscoelastic  effects  should  also  be  taken 
into  account.  Only  elastic  response  is  considered 
here.  Then  (27)  implies  that  pu(p)  differs  from 
Ph(p)  by  a  term  of  order  (p  —  l)3  unless  Co  =  0, 
which  is  generally  not  the  case.  If  cq  and  Ph(p) 
are  known,  then  (27)  provides  an  approximation 
to  pu  to  within  an  error  of  order  (p—  l)3.  The  rel¬ 
ative  error  in  approximating  pu  by  p*  is  of  order 
p  —  1  and  can  be  estimated  by  using  (28). 

In  a  pressure/shear  (or  oblique  plate  impact) 
test,  a  longitudinal  wave  propagates  into  the  un- 
deformed  material  bringing  it  to  a  state  of  uni¬ 
axial  strain,  and  a  slower  shear  wave  propagates 
into  this  uniaxial ly  strained  material.  These  tests 
yield  both  tfi(p)  and  the  shear  wave  speed  Ustu 
(and  hence  Gu)  as  a  function  of  p  or  &\.  If  the 
shear  wave  travels  at  the  acceleration  wave  speed 
then  (17),  (19),  and  (21)  provide  exact  formulas 
for  the  shear  stress  r,  the  pressure  and  the 
bulk  modulus  Ku  in  uniaxial  strain  as  a  func¬ 
tion  of  p  or  <ti  .  These  formulas  appear  to  have 
gone  unnoticed,  however.  Instead,  it  is  usually 
assumed  that  Ku  zs  Lu  —  |GU.  This  approxi¬ 
mate  relation,  together  with  (20),  is  then  inte¬ 
grated  to  give  pu  as  a  function  of  p\  cf,  Gupta 
(4,6),  Conner  (5),  and  Aidun  k  Gupta  (7)+  For 
fused  silica  in  the  strain  range  0  <  p  -  1  <  0.076, 
the  response  is  elastic  and  the  shear  wave  speed 
decreases  with  p;  cf.  Conner  (5).  In  this  strain 
range  the  shear  wave  is  an  acceleration  wave  (cf. 
also  Abou-Sayed  k  Clifton  (8))p  so  we  may  ap¬ 
ply  the  results  of  §4.  Using  (21)  and  Conner's 
data,  we  find  that  at  a  strain  of  p  —  1  a  0.076  the 
estimate  Ku  ss  Lu  —  is  low  by  about  29%. 

Whether  the  shear  wave  in  a  pr ess u re/shear 
test  is  an  acceleration  wave  or  a  shock  wave  de¬ 
pends  on  the  nonlinear  elastic  response  of  the 
material;  cf.  Davison  (9).  The  shear  modulus 
Gu  in  §4  is  defined  in  terms  of  the  accelera¬ 
tion  wave  speed  or  equivalently,  in  terms 

of  the  speed  of  a  plane  infinitesimal  sinusoidal 
shear  wave;  cf.  |2,  If  a  shear  shock  with  speed 
U  can  propagate  in  the  uniaxially  strained  mate¬ 
rial  and  if  we  set  G  M  pU2,  then  the  formulas  in 
§4  hold  approximately  when  Gu  is  replaced  with 


G,  Also  note  that  if  U  >  UsiU  (as  standard  stabil¬ 
ity  arguments  would  imply),  then  0  >  GUi  and 
(17)  and  (19)  imply  that  r  <  ^(p3  —  1)G  and 
Pu  >  —  §(p2  “  1)<5  in  compression  (p  >  1), 

We  conclude  with  a  brief  discussion  of  thermo¬ 
dynamic  effects,  which  have  been  neglected  up 
to  this  point.  If  a  thermoelastic  material  con¬ 
ducts  heat  by  Fourier's  law  [resp.  is  a  nonconduc¬ 
tor],  then  a  longitudinal  acceleration  wave  prop¬ 
agates  at  the  isothermal  [resp.  adiabatic]  wave 
speed.  However,  the  formula  (15)  for  the  speed 
of  a  shear  acceleration  continues  to  hold  in  ei¬ 
ther  case;  cf.  Bowen  k  Wang  (ID).  In  fact,  it  can 
be  shown  that  (15)  holds  even  if  heat  conduc¬ 
tion  is  governed  by  Cattaneo's  equation,  which 
prohibits  instantaneeous  propagation  of  thermal 
disturbances.  Thus  the  formulas  (17)  and  (19)  for 
the  shear  stress  and  the  pressure  continue  to  hold. 
In  particular,  they  are  valid  when  the  state  of  uni¬ 
axial  strain  has  been  achieved  by  shock  loading. 

REFERENCES 

1 .  Stheidier,  M.,  “On  the  coupling  of  pressure  and 
deviatoric  stress  in  hyperelastic  materials,'1  in  Pro¬ 
ceedings  of  the  13th  Army  Symposium  on  Solid  Me¬ 
chanics ,  1993,  pp.  539-550. 

2*  Truesdell,  C„  and  Noll,  W„  The  Non- Linear  Field 
Theories  of  Mechanics,  Berlin:  Springer,  1905, 
S48,  573-75* 

3.  Wang,  C.-C.,  and  Thiesdell,  C<,  Introduction  to 
Rational  Elasticity,  Leyden,  Noordhoff,  1973, 
$VI.5* 

4.  Gupta,  Y,  M.,  A  Appl.  Phys 51,  5352-5361 
(1980). 

5.  Conner,  M.  P.,  “Shear  Wave  Measurements  to  De¬ 
termine  the  Nonlinear  Elastic  Response  of  Fused 
Silica  Under  Shock  Loading,"  M.  S.  thesis,  Wash¬ 
ington  State  University,  Pullman,  1988. 

6.  Gupta,  Y,  M.,  A  Geophys.  Res.,  88,  4304^312 
(1983). 

7.  Aidun,  J,  B.,  and  Gupta,  Y,  M.,  Geophys .  Res. 
Utters,  16,  191-194  (1989). 

8.  Abou-Sayed,  A.  S,,  and  Clifton,  R.  J.,  J.  Appl. 
Phys.,  47,  1762-1770  (1976). 

9.  Davison,  L.,  A  Meeh*  Phys.  Solids,  14,  249-270 
(1966)* 

10.  Bowen,  R.  M.,  and  Wang,  C.-C.,  Arch .  Anlinna/ 
Mech.  Anal.,  41,  287-318  (1971). 


478 


CP429,  Shock  Compression  of  Condensed  Matter- 1997 
edited  by  Schmidt/Dandekar/Forbes 

1998  The  American  Institute  of  Physics  l-56396-738-3/98/$15.00 


RESPONSE  OF  NONLINEAR  ELASTIC 
SOLIDS  TO  OBLIQUE  PLATE  IMPACT 

Mike  Scheidler 

U.S.  Army  Research  Laboratory,  APG,  Maryland  21005-5066 

We  give  a  theoretical  analysis  of  the  nonlinear  elastic  response  in  the  interior  of  the  target  in  an 
oblique  plate  impact  test,  prior  to  interactions  with  any  reflected  waves.  Approximate  relations  are 
derived  between  the  changes  in  stress,  strain,  particle  velocity,  and  wave  speeds  across  the  shear  wave 
by  neglecting  shear  strain  terms  of  order  six  in  the  internal  energy  function.  The  results  are  valid 
for  any  finite  (elastic)  longitudinal  strain  ahead  of  the  shear  wave.  They  apply  to  isotropic  materials 
and  to  appropriately  aligned  transversely  isotropic  and  orthotropic  materials. 


INTRODUCTION 


PRELIMINARIES 


In  an  oblique  plate  impact  test  a  flyer  plate, 
inclined  relative  to  the  axis  of  the  projectile,  im¬ 
pacts  a  parallel  target  plate,  which  we  assume  is 
at  rest  and  stress  free.  Let  (X,  Y,  Z)  denote  the 
Cartesian  coordinates  of  a  material  point  in  this 
natural  reference  state,  with  the  X-axis  normal 
to  the  target  face.  The  particle  velocity  vF  im¬ 
parted  to  the  face  of  the  target  has  nonzero  com¬ 
ponents  both  normal  and  parallel  to  the  target 
face.  The  coordinate  axes  are  oriented  so  that 
vF  =  (?/F,uF,0),  with  the  normal  (X)  component 
uF  >  0  and  the  transverse  (Y)  component  vf  >  0. 
The  impact  is  assumed  weak  enough  that  the  tar¬ 
get  response  is  elastic.  We  assume  the  target  is 
homogeneous,  and  restrict  attention  to  points  in 
the  interior  and  times  prior  to  the  arrival  of  waves 
reflected  from  free  surfaces  or  material  interfaces. 
Then  the  motion  is  independent  of  the  Y  and  Z 
coordinates.  However,  for  a  general  anisotropic 
material,  the  motion  may  have  nonzero  compo¬ 
nents  in  the  Z-direction  at  points  not  on  the  im¬ 
pact  face.  For  isotropic  materials  as  well  as  some 
anisotropic  materials  (appropriately  aligned  with 
the  coordinate  axes),  there  is  no  motion  in  the 
Z-direction.  This  is  the  case  considered  here. 


If  (x,y,z)  is  the  position  at  time  t  of  the 
material  point  initially  at  (X,  Y,  Z),  then  x  = 
X  +  di(X,f),  y  =  Y  +  d2(X,t),  *  =  Z.  The 
normal  and  transverse  components  of  the  par¬ 
ticle  velocity  are  u  —  dx/dt  and  v  =  dy/dt. 
Let  P  and  T  denote  the  1st  Piola-Kirchhoff  and 
Cauchy  stress  tensors.  Then  a  =  -Pxx  =  —Txx 
and  e  =  1  —  dx/dX  are  the  normal  or  longitu¬ 
dinal  components  of  stress  and  strain,  taken  pos¬ 
itive  in  compression;  and  r  =  —Pyx  —  —Tyx 
and  7  =  -dy/dX  are  the  shear  stress  and  shear 
strain.  The  signs  have  been  chosen  so  that  e,  7, 
cr,  r,  as  well  as  u  and  v,  should  be  nonnegative 
for  the  impact  problem  considered  here.  Since  e 
and  7  are  the  only  nonzero  strain  components, 
the  internal  energy  e  per  unit  mass  depends  only 
on  e,  7,  and  the  entropy  per  unit  mass  s.  Then 

-v  x  de  de 

e  =  e(e,7, s) ,  cr  =  p0— ,  r  =  p0—,  (1) 

where  p0  is  the  density  in  the  natural  state.  Heat 
conduction  is  neglected  in  the  analysis,  so  that 
response  is  isentropic  except  for  entropy  jumps 
across  shocks.  The  equations  of  motion  are 


du  da  dv  dr 

Po~dt+dX~  ’  Po~dt  +  dX=°' 


(2) 
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The  strain  rate  and  velocity  gradient  satisfy 
du 


de 

m  +  ax  ~  ’ 


dj  dv 

~dt+dx=:  ' 


(3) 


Smooth  solutions  are  governed  by  (2)  and  (3), 
with  a  and  r  given  by  (1)  and  s  constant;  this  is 
a  quasilinear  system  in  the  four  unknowns  u,  v, 
£,7.  The  (Lagrangian)  characteristic  wavespeeds 
are  the  roots  ±U  and  ±V  of 


2p0 


_  da  dr  If  da  <9r\2  7 <9r\2 

+  ^7  y  V  9e  <9y 7  /  ’  (4) 


where  the  fast  wave  speeds  ±1/  correspond  to 
the  +  sign  and  the  slow  wave  speeds  ±V  to  the 
-  sign  in  (4).  The  constitutive  inequalities 


da  dr  da  dr  /<9r\2 

de  >  d'y  >  ’  de  dj  >  \de ) 


(5) 


guarantee  that  there  are  four  real  distinct  charac¬ 
teristic  wave  speeds  (i.e.,  U  >  V  >0),  and  that 
the  fast  waves  are  primarily  longitudinal.  Waves 
with  negative  speed  would  be  generated  by  re¬ 
flections  and  hence  are  not  considered  here.  The 
right  inequality  in  (5)  is  equivalent  to 


d2e  d2e  (  d2e  \2 
de2  dj2  >  \  ded^f ) 


(6) 


which  implies  that  e  is  a  strictly  convex  function 
of  e  and  7. 

The  one-dimensional  plane  waves  for  this 
oblique  impact  problem  are  centered  simple  waves 
and/or  centered  shocks.  For  nonlinear  isotropic 
elastic  solids,  some  of  the  earlier  papers  on  this 
problem  are  Bland  (1),  Davison  (2),  and  Abou- 
Sayed  &  Clifton  (3).  Isotropy  implies  that  e  is  an 
even  function  of  7  for  fixed  e  and  s,  and  that  the 
fast  wave  is  purely  longitudinal,  i.e.,  there  are  no 
changes  in  v,  7,  or  r.  The  fast  wave  brings  the 
material  to  an  intermediate  state  of  uniform  uni¬ 
axial  strain.  The  slow  wave  propagates  into  this 
uniaxially  strained  material.  The  slow  wave  is 
primarily  a  transverse  or  shear  wave.  However, 
nonlinear  elastic  effects  cause  2nd  order  changes 
in  u,  e,  and  a  across  the  shear  wave.  The  case 
where  both  waves  are  simple  waves  is  illustrated 
in  Fig.  1. 


FIGURE  1.  Simple  wave  solutions  for  oblique  plate  impact. 
Quanties  in  th4e  intermediate  state  (ahead  of  the  shear  wave) 
and  final  state  (behind  the  shear  wave)  were  denoted  by  I  and  F 
subscriptes,  respectively. 

For  anisotropic  materials  there  may  be  motion 
in  the  Z  direction.  Even  if  there  is  no  motion 
in  the  Z  direction,  the  fast  wave  need  not  be 
purely  longitudinal.  However,  if  the  material  is 
orthotropic  relative  to  the  X,  F,  Z  axes,  or  trans¬ 
versely  isotropic  about  the  X  or  Z  axis,  then  it 
can  be  shown  that  there  is  no  motion  in  the  Z  di¬ 
rection,  and  that  e  is  an  even  function  of  7.  Fur¬ 
thermore,  it  can  be  shown  that  this  last  property 
implies  all  of  the  results  discussed  previously  for 
the  isotropic  case.  With  these  restrictions,  most 
of  the  analysis  for  isotropic  materials  in  references 
(l)-(3)  can  be  extended  to  the  anisotropic  case. 
Since  e  is  an  even  function  of  7  for  fixed  e  and  5, 

e  «  f{e,  s)  +  \  g{e,  s)  72  +  £  h{e,  s)  j4.  (7) 

Here  and  below,  a  «  b  means  a  =  b  to  within  an 
error  of  order  7n.  Note  that  no  approximations 
are  made  for  the  dependence  of  e  on  e  or  s.  In  par¬ 
ticular,  the  results  in  this  paper  are  valid  for  ar¬ 
bitrary  finite  (elastic)  longitudinal  strains  ahead 
of  the  shear  wave.  We  study  only  the  shear  wave 
in  this  paper,  and  consider  only  the  case  where  it 
is  a  simple  wave  (as  opposed  to  a  shock).  s0  and 
sl  denote  the  entropy  in  the  natural  and  interme¬ 
diate  states  .  If  the  longitudinal  wave  is  simple, 
Sj  =  s0 ;  if  it  is  a  shock,  s1>s0. 
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SIMPLE  SHEAR  WAVES 

Set  f(e)  -  /(e, s,)>  g(e)  =  g{e, s,),  and  h(e)  = 
h(e,  Sj).  Since  the  deformation  in  the  simple  shear 
wave  is  isentropic,  (7),  (1),  and  (4)  imply 


e  w /(e) +  |  0(e)  ^  h(e)  ^ ,  (8) 

y  =  ||  «/'(£)  +  Iff' (e)72,  (9) 

—  =  «  5(e)  7+  £M073>  (10) 

Po 


where  a  prime  denotes  differentiation  with  respect 
to  e.  Note  that  o  and  V  are  even  functions  of  7 
for  fixed  e,  while  r  is  an  odd  function  of  7. 

In  the  intermediate  uniaxially  strained  state 
ahead  of  the  shear  wave,  7  =  r  =  v  =  0,  and  by 
(9)— (1 1) ,  (4),  and  (5)i,  a1  =  />„/'(£,)  and 

7  =  /"(£l)  >  g(£l)  =  V,2.  (12) 


The  rate  of  change  of  the  speed  V1  of  the  shear 
wave  front  (which  is  a  transverse  acceleration 
wave)  with  respect  to  the  uniaxial  strain  ahead 
of  the  wave  is 


r,_rfv,  _  g'(Q 
de1  2  \AKeJ 


(13) 


In  the  centered  simple  shear  wave,  F  —  X/t\ 
and  £,  a,  u,  7,  r  &  v  are  constant  along  the 
straight-line  characteristics  (i.e.,  functions  of  X/t 
only).  The  characteristic  form  of  the  quasilinear 
system  yields  the  differential  relations 


du 

dv 


do- 

dr 


de  _ 
d^y  <9a 
de 


dr 

de 


-P0V2 


dr 

dj 


-pav\ 


(14) 

(15) 


7,  r  &  v  are  strictly  increasing  with  passage  of 
the  wave,  and  any  two  of  these  variables  is  an  odd 
function  of  the  other.  £,  <7,  u  &  V  are  even  func¬ 
tions  of  any  one  of  the  variables  7,  r  or  v\  hence 
the  total  derivatives  in  (14)  are  0  at  the  wavefront, 


where  7  =  r  =  v  =  0.  Indeed,  do /de  —  p0V2  >  0 
by  (5)i,  while  by  (10),  dr/de  =  0  at  the  wave- 
front.  If  dr/de  ^  0  for  7^0,  then  (15)  holds 
with  v,7 hr  replaced  by  u,eho,  respectively. 
Inequality  (5)3  places  no  restriction  on  the  sign  of 
dr/de ;  if  it  changes  sign  within  the  shear  wave, 
then  by  (14),  u,  o  &  e  do  not  vary  monotonically 
through  the  wave. 

Next,  we  quantify  the  above  statements  by  de¬ 
riving  approximate  relations  which  hold  through¬ 
out  the  shear  wave.  On  using  (9)-(13)  in  (14)3 
and  noting  that  e  =  el  when  7  =  0,  we  obtain 

£Ksl-\Q{el)  72,  (16) 


5'  2  V,  V/ 


~  f"  -  ~  tt2  1/2 


(17) 


9  -  u2  -v; 

On  substituting  (16)  into  (11)  and  using  (12), 
(13),  (17),  and  V  =  V1  when  7  =  0,  we  obtain 

V 


a~VI2-R(eI)72,  V^V,-^^2,  (18) 


2V, 


2R(6i)  =  -MO 


12  V2 


(19) 


(V,T  -  MO- 


U2  -  V2 

From  (15)i,  (18)2,  and  v  =  0  when  7  =  0,  we  get 


5  R(e,)  3 

u  ~  V,  7  QyT-  73, 


which  may  be  inverted  to  give 

5  v  R(e  1)  (  v 

+  'elf  \vi 

Then  (18)2,  (21),  and  (15)2  imply 

vAv,  R{c’)l'v 


T  5 


Vtv- 


2  V,  VV, 
R(e,) 


V  . 


Pc  ‘  6  V, 

Similar  arguments  yield  the  relations 

£«£,  -  |Q(e,)  (u/V,)2, 
«  cr,  -  1  <5(e,)p0n2, 
u  «  u,  -  |  Q(e,)  v2 /Vl. 


(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
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Dropping  the  v3  terms  in  (21)  and  (23)  and  elim¬ 
inating  Q(ez)  from  (24)-(26)  yields 

7  w  v/Vj ,  e-eI  «  (u-uJ/Vj,  (27) 

3  4 

t  «  PoV.v,  a  -  al  «  p0VM  -  u,) .  (28) 

More  accurate  material-independent  approxima¬ 
tions  for  7  and  r  can  be  obtained  by  solving  (22) 
for  Rfa)  and  substituting  into  (21)  and  (23): 

7  «  [1  +  WJVr-  1  )(v/t02]v/V;,  (29) 

(30) 

Since  V  must  decrease  with  passage  of  the 
wave,  (18)  implies  that  Rfa)  >  0  is  necessary 
for  the  centered  shear  wave  to  be  simple  (as  op¬ 
posed  to  a  shock)  for  sufficiently  weak  impacts, 
while  R(e1)>  0  is  sufficient.  It  can  be  shown  that 
R(e J  <  0  is  sufficient  for  a  shear  shock.  There 
are  no  restrictions  on  Q(er),  however.  By  (17) 
and  (12),  we  see  that  Q{el)  and  VT'  have  the  same 
sign.  Hence,  (24)-(26)  imply  that  if  Q(ex)  or  V/ 
is  negative  [resp.  positive],  then  £,cr, u  increase 
[resp.  decrease]  across  the  shear  wave  for  suffi¬ 
ciently  weak  impacts;  if  Q{e^)  or  V/  is  zero,  then 
£,  o',  u  are  constant  to  within  an  error  of  order  74. 

Solving  (18)2  for  7  in  terms  of  V  shows  that  7 
has  a  square  root  singularity  at  V  =  V1,  i.e.,  at  the 
wavefront;  likewise  so  do  r  and  v.  Since  V  =  X/t, 
the  strain  rate,  stress  rate,  and  particle  acceler¬ 
ation  are  infinite  at  the  wavefront,  as  noted  by 
Abou-Sayed  &  Clifton  (3)  for  isotropic  materials. 

At  the  rear  of  the  shear  wave  (i.e.,  at  the  final 
state),  (29)  and  (30)  reduce  to 


7f 


-  +  ~  !!i4(i  +  ±a)  — 

3  +  3vJv;  1  3^V,’ 


(31) 


rF  «  Po(f^i  +  IHK  «  (1  -  |A s)p0VlVf,  (32) 


where  As  =  VJVF  -  1  =  V1tR(X)/X,  and 
tn(X)  —  X/Vf -  X/Vl  is  the  rise  time  of  the  shear 
wave  at  X.  By  (22), 


K 


Vn 


1  4  Rte.) 
'v^Tvf  \X 

6(7)2  MO)' 


U‘ 


2VnJ 


Vn 


2 


(33) 


The  top  formula  in  (33)  includes  dependence 
of  the  rise  time  on  the  uniaxial  strain  el  ahead  of 
the  shear  wave.  The  bottom  formula  neglects  this 
dependence;  it  is  obtained  by  setting  el  =  0  and 
using  (19),  and  contains  errors  of  order  7F4  and 
£j7f2.  Here  U0  and  V0  are  the  longitudinal  and 
shear  wave  speeds  in  the  natural  state,  and  V0f 
denotes  the  value  of  VT'  at  el  =  0.  Thus  U0&tV0, 
Vq,  and  h( 0)  are  2nd,  3rd,  and  4th  order  elastic 
constants,  respectively.  An  analogous  formula  in 
Abou-Sayed  &  Clifton  (3),  eqn.  (46),  is  off  by  a 
factor  of  3/2  (the  result  of  neglecting  changes  in  e 
across  the  shear  wave),  and  missing  the  h( 0)  term 
(the  result  of  neglecting  4th  order  strain  terms  in 
the  internal  energy). 

From  (22)  and  (19)  we  obtain  an  approxima¬ 
tion  for  the  4th  order  coefficient  h  in  (8)  evaluated 
at  the  uniaxial  strain  e1  ahead  of  the  shear  wave: 


Me  1) 


12  V,2 

u?  -  7 


(7)  -  47 


K 


(34) 


V/  can  be  estimated  by  differentiating  the  V1  (eI ) 
curve  obtained  from  a  series  of  oblique  plate  im¬ 
pact  tests,  or  by  solving  (26)  and  (17): 


V* 


)  (7  -  7)-  (35) 


When  this  is  substituted  into  (34),  we  obtain 


h{e. 


,)  »  12  (^70  ^  ~  V'!>  *7 


(36) 


The  right-hand  sides  of  (35)  and  (36)  involve  par¬ 
ticle  velocities  and  wave  speeds  only.  These  could 
be  measured  in  a  single  test,  at  least  if  the  longi¬ 
tudinal  wave  is  simple.  If  the  longitudinal  wave 
is  a  shock  with  speed  U,  then  C/r2  can  be  approx¬ 
imated  by  U?  «  2U2  -  U0 2  to  within  an  error  of 
order  £j2,  or  calculated  from  (12)i  if  /  is  known. 
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UNIVERSAL  RELATIONS  FOR  PRESSURE-SHEAR 
WAVES  IN  NONLINEAR  ELASTIC  SOLIDS 

Mike  Scheidler 

U.S.  Army  Research  Laboratory ,  APG,  Maryland  21005-5069 

The  centered  shock  waves  or  simple  waves  generated  by  oblique  plate  impact  of  nonlinear  elastic 
solids  are  studied.  Attention  is  restricted  to  materials  for  which  the  fast  wave  is  purely  longi¬ 
tudinal.  This  includes  isotropic  materials  and  appropriately  aligned  orthotropic  or  transversely 
isotropic  materials.  The  emphasis  is  on  the  derivation  of  exact  or  approximate  relations  between 
stress,  strain,  particle  velocities,  and  wave  speeds  which  are  universal  in  the  sense  that  they  do 
not  involve  material  constants  or  material  functions.  For  isotropic  materials,  universal  relations 
for  the  difference  between  the  longitudinal  and  lateral  stresses  are  derived. 


INTRODUCTION 

Consider  a  homogeneous  target  plate  initially 
at  rest  and  stress-free.  (X,Y,Z)  are  the  Carte¬ 
sian  coordinates  of  a  typical  point  in  this  natural 
reference  state,  with  the  X-axis  normal  to  the 
target  face.  In  an  oblique  plate  impact  test,  the 
particle  velocity  (wF,t;F,0)  imparted  to  the  face 
of  the  target  by  the  flyer  plate  has  a  normal  (X) 
component  uF  >  0  and  a  transverse  (Y)  compo¬ 
nent  vf  >  0.  The  impact  is  assumed  weak  enough 
that  the  response  of  the  flyer  and  target  is  elas¬ 
tic.  We  restrict  attention  to  points  in  the  target 
interior  and  to  times  prior  to  the  arrival  of  waves 
reflected  from  free  surfaces  or  material  interfaces. 
Then  the  motion  consists  of  plane  pressure- shear 
waves  propagating  in  the  X  direction.  For  nonlin¬ 
ear  isotropic  elastic  solids,  this  problem  has  been 
studied  by  Bland  (1),  Davison  (2),  and  Abou- 
Sayed  &  Clifton  (3).  As  noted  by  Scheidler  (4), 
many  of  their  results  remain  valid  if  the  material 
is  orthotropic  relative  to  the  X,  Y,  Z  axes,  or  if  it 
is  transversely  isotropic  about  either  the  X  axis 
or  the  Z  axis. 

Assuming  one  of  the  above  material  symmetry 
conditions,  the  motion  has  the  general  form 

x=X+d(X,t),  y  =  Y+d(X,t),  z  =  Z,  (1) 


where  (x,y,z)  is  the  position  at  time  t  of  the  ma¬ 
terial  point  initially  at  (XjY,Z).  Let  P  and  T 
denote  the  1st  Piola-Kirchhoff  and  Cauchy  stress 
tensors.  Then  the  normal  or  longitudinal  com¬ 
ponents  of  stress  and  strain,  taken  positive  in 
compression,  are 

dx 

°x  =  ~Txx  =  ~Pxx,  (2) 


The  lateral  stress,  taken  positive  in  compression, 
is 


<JY  —  Tyy 


~^YY  T  Lx  Y 

1  -  £ 


(3) 


The  shear  stress  and  shear  strain  are 


r  =—T  —-T  =-P  'Y  = 

1 XY  —  Xxy  *YX  x  YX  5  I  - 


dy 

dx 


(4) 


The  signs  have  been  chosen  so  that  e,  7,  ax ,  aY , 
and  rxr,  as  well  as  the  normal  and  transverse 
components  of  the  particle  velocity,  u  —  dx/dt 
and  v  =  dy/dt ,  should  be  nonnegative  for  the 
impact  problem  considered  here. 

Since  e  and  7  are  the  only  nonzero  strain  com¬ 
ponents,  e  =  e(e,7,s),  where  e  and  s  are  the  in¬ 
ternal  energy  and  entropy  per  unit  mass.  Heat 
conduction  is  neglected  in  the  analysis,  so  the 
response  is  isentropic  except  for  entropy  jumps 
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across  shocks.  p0  denotes  the  density  in  the  nat¬ 
ural  state.  The  material  symmetry  restrictions 
imply  that  e  and  ax  =  p0de/de  are  even  functions 
of  7  for  fixed  e  and  s,  whereas  rXY  =  p0de/d 7  is 
an  odd  function  of  7.  These  even/odd  properties 
result  in  qualitative  differences  in  the  behavior  of 
longitudinal  and  shear  waves. 

Two  centered  waves  are  generated  on  impact. 
The  material  symmetry  restrictions  imply  that 
the  fast  wave  is  purely  longitudinal,  and  so  brings 
the  material  to  an  intermediate  state  of  uniform 
uniaxial  strain  into  which  the  slower  shear  wave 
propagates.  The  shear  wave  brings  the  mate¬ 
rial  to  a  final  state  of  uniform  compression  and 
shear.  Quantities  in  the  natural,  intermediate, 
and  final  states  are  denoted  by  0,  I,  and  F  sub¬ 
scripts,  respectively.  To  simplify  the  notation, 
we  let  a  =  ax  and  r  =  rXY.  In  particular, 
7j  =  Tj  =  vl  =0.  The  shear  wave  is  often 
called  a  quasi-transverse  wave  since  it  also  in¬ 
duces  changes  in  the  longitudinal  quantities  e, 
(j,  and  u.  For  the  classes  of  materials  considered 
here,  these  changes  are  2nd  order  in  ,  whereas 
rF  and  vF  are  1st  order  in  yF . 

For  sufficiently  small  strains,  the  centered 
waves  are  either  shocks  or  simple  waves,  depend¬ 
ing  on  the  nonlinear  elastic  properties  of  the  ma¬ 
terial.  The  case  where  both  waves  are  simple  is 
illustrated  in  Fig.  1. 


FIGURE  1.  Simple  waves  generated  by  oblique  plate 
impact.  Quantities  in  the  intermediate  state  (ahead  of 
the  shear  wave)  and  final  state  (behind  the  shear  wave) 
are  denoted  by  I  and  F  subscripts,  respectively. 


The  Lagrangian  acoustic  (or  characteristic  or 
acceleration)  longitudinal  and  shear  wave  speeds 
are  denoted  by  U  and  V,  respectively.  They  are 
even  functions  of  7  for  fixed  e  and  s.  For  the 
case  of  a  simple  shear  wave,  this  implies  that  the 
time  histories  of  7,  r,  and  v  have  square  root 
singularities  at  the  shear  wave  front;  cf.  Abou- 
Sayed  &  Clifton  (3)  for  the  isotropic  case  and 
Scheidler  (4)  for  the  more  general  case  consid¬ 
ered  here.  On  the  other  hand,  for  a  simple  lon¬ 
gitudinal  wave  the  time  histories  of  e,  a,  and 
u  are  smooth  (and  hence  approximately  linear) 
throughout  the  wave,  and  such  waves  are  often 
referred  to  as  ramp  waves.  Since  the  analysis  of 
the  longitudinal  wave  is  standard,  the  remainder 
of  this  paper  focuses  on  the  shear  wave.  The 
results  are  valid  for  arbitrary  finite  (elastic)  lon¬ 
gitudinal  strains  sl  ahead  of  the  shear  wave.  Due 
to  limitations  of  space  some  proofs  are  omitted. 
For  the  case  of  a  simple  shear  wave,  derivations  of 
some  of  the  results  can  be  found  in  Scheidler  (4). 

SHEAR  SHOCKS 

In  this  section  we  consider  the  case  where  the 
shear  wave  is  a  shock.  The  jump  conditions  are 


II 

ur  —  u T 

£  £  —  1 

F  1  v  ’ 

(5) 

TF  =  PoV  > 

ov  -  -  «i) , 

(6) 

where  ¥  is  the  Lagrangian  shear  shock  speed. 
These  relations  imply 

Tf  =  pfy2%  ,  <TF  -  <7,  =  p0V2 (eF  -  e,).  (7) 

¥,  VF ,  and  Vl  all  differ  by  terms  of  order  7F2 .  Sta¬ 
bility  arguments  imply  that  VF  >  ¥  >  Vl,  i.e.,  the 
shock  wave  speed  is  bounded  below  by  the  acous¬ 
tic  wave  speed  ahead  of  the  shock  and  above  by 
the  acoustic  wave  speed  behind  the  shock.  ¥  may 
be  approximated  by  the  average  of  these  acoustic 
wave  speeds:  ¥  =  ^V1  -f  \VF  4-  0( 7F2).  Such  an 
approximation  will  be  abbreviated  by 

y*kv1  +  \vr,  (8) 

where,  in  general,  a  £  b  means  that  a  —  b  to 
within  an  error  of  order  7Fn. 
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An  interesting  result,  which  has  no  analog  for 
longitudinal  shocks,  is  that  the  approximation  (8) 
may  be  improved  simply  by  changing  the  weight¬ 
ing  factors: 


(9) 


This  is  a  consequence  of  the  fact  that  r  is  an  odd 
function  of  7.  (9)  is  equivalent  to  either  of  the 
relations 

(10) 

H-v  «§(VF-V0.  (11) 

Bland(1965,  p.764)  obtained  (10)  for  the  special 
case  of  isotropic  materials  with  infinitesimal  uni¬ 
axial  strain  ahead  of  the  shear  wave.  Using  the 
fact  the  sum  of  the  weighting  factors  in  (9)  is  1, 
it  may  be  shown  that  (9)— (11)  are  equivalent  to 
analogous  formulas  for  the  squares  of  the  wave 
speeds: 


ly  2  ,  ly  2 

31  '  3  F  J 

(12) 

tG^-v,2), 

(13) 

f  (Vf2  -  V,2) . 

(14) 

yF2-^ 


SIMPLE  SHEAR  WAVES 

Although  s,  7,  etc.  vary  continuously  through¬ 
out  a  centered  simple  wave,  it  is  well-known 
that  the  net  changes  across  such  a  wave  may  be 
approximated  by  the  jump  conditions  across  a 
shock,  with  the  shock  speed  replaced  by  the  speed 
of  the  simple  wave.  In  particular,  for  a  centered 
simple  shear  wave  we  have  the  following  approx¬ 
imate  analogs  of  the  jump  conditions  (5)— (7): 


In  the  above  approximations  the  speed  V1  of  the 
simple  wave  front  may  be  replaced  by  the  speed 
VF  of  rear  of  the  wave.  More  accurate  approxi¬ 
mations  can  be  obtained  by  replacing  VT  in  (15)- 
(17)  by  appropriate  weighted  averages  of  Vl  and 
VL .  For  the  transverse  variables  we  have 


5 

7f  ~  2 


V  4-  ~V 

VI  ^  SVF 


(18) 


Tp  +  KK 

«  P0(|^I  +  |^f)27f 

«Po(|^2  +  ^f2)7f- 
For  the  longitudinal  variables  we  have 
6  uF  —  ux 

e*~El~  f^  +  IVp  5 

<JF  —  (Jj  ~  P0  (  2  ^1  f  2^f)  (^f  ~  ^i) 
~  ^0(2^1  +  2^f)2  (£f  “  £i) 

«Po(^2  +  K2)(eF-^) 


(19) 


(20) 


(21) 


APPROXIMATIONS  FOR  9VJdel 

Since  V  =  in  the  intermediate  uni- 

axially  strained  state  where  7  =  =  0  it  follows 

that  the  acoustic  shear  wave  speed  V  depends 
only  on  the  uniaxial  strain  e  =  e1  and  entropy 
s  =  5i;  V[  =  V’(eI,0,5I).  ThendVjdSj  is  the  rate 
of  change  of  the  acoustic  shear  wave  speed  with 
respect  to  uniaxial  strain,  evaluated  at  the  inter¬ 
mediate  uniaxially  strained  state.  If  the  oblique 
plate  impact  generates  a  shear  shock,  then 


2  aF 


uF  —u. 


dVr 


U 2  -  ¥2  8eI  ’ 


(22) 


where  the  top  formulas  follow  immediately  from 
(5)-(6) ,  and  the  bottom  formula  can  be  obtained 
by  expanding  the  stress  quotient  in  powers  of  7F . 


4  Uf  —  U, 

p  p  r  1 

(15) 

Similarly,  if  tl 

v 

fcF  ~  y  ? 

eF  -e, 

Po^I  ^F  5 

0-F  -  (Tj  W  pnVx  (uf  - 

«i)»  (16) 

7f2 

O'f  -  «  PoVi2(£r  - 

-  e,).  (17) 

P»V? 


£  y  —  «, 

>  ~  vi  ..2 


-K 


dK 


U‘2  -  V*  ds , 


(23) 


where  the  top  formulas  follow  from  (15)-(16).  In 
(23)  the  speed  Vt  of  the  simple  wave  front  may 
be  replaced  by  the  speed  V¥  of  rear  of  the  wave. 
If  we  set 

{¥ ,  for  shear  shock 

V1  or  UF  ,  for  simple  shear  wave  (^4) 

and  solve  (22)  or  (23)  for  dVjde^  we  obtain 
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Up  «  -(U,2  -  v2)(uf  -  u^/v? 

(26) 

*-P,V{U?-V2){<Tf-ai)lT?. 

Simultaneous  measurement  of  longitudinal 
and  transverse  particle  velocities  and  wave  speeds 
in  oblique  plate  impact  tests  have  been  per¬ 
formed  by  Gupta  and  coworkers  using  in-material 
electromagnetic  particle  velocity  gauges;  cf.  (5). 
Such  measurements  could  be  used  to  estimate 
dV1/ds1  by  means  of  the  relation  (25)  i. 

ISOTROPIC  ELASTIC  SOLIDS 

We  return  to  the  original  notation  of  crx  and 
txy  for  the  normal  and  shear  stresses.  Since  the 
lateral  stress  aY  does  not  appear  in  the  momen¬ 
tum  balance  equations,  the  wave  analysis  yields 
no  information  on  aY ,  at  least  for  the  anisotropic 
solids  considered  up  to  this  point.  However,  isot¬ 
ropy  places  restrictions  on  aY ,  as  we  now  show. 

F  =  dx/dX  is  the  deformation  gradient.  For 
an  isotropic  elastic  solid,  the  Cauchy  stress  ten¬ 
sor  T  commutes  with  B  —  FFr :  TB  —  BT. 
For  the  pressure-shear  motion  (1),  the  only  non¬ 
trivial  restriction  that  follows  from  this  result  is 

Txy  (Bxx  ~Byy)  =  Bxy  (Txx  -  Tyy  )•  For  7  =  0 
(i.e.,  uniaxial  strain)  this  reduces  to  0  =  0,  but 
for  7  ^  0  it  implies  the  universal  relations 


ference  between  the  normal  and  lateral  stresses  in 
the  intermediate  uniaxially  strained  state  ahead 
of  the  shear  shock: 


2s- 


PoV 2 


2 

-1 


Pv\ 

(28) 


For  a  simple  shear  wave,  a  similar  analysis  using 
(15)i~(17)i  also  yields  (28).  Here  and  below,  we 
drop  the  I  subscript  since  all  quantities  are  eval¬ 
uated  in  the  intermediate  state.  (28)2  was  first 
derived  by  Scheidler  (6)  by  a  different  method; 
V  =  (1  —e)V  is  the  Eulerian  acoustic  shear  wave 
speed  and  p  the  density.  (28)  also  yields  formulas 
for  the  pressure  p  and  maximum  shear  stress  rmax 
in  the  uniaxially  strained  state,  since 


2  Tmax  &x  ?  P  3  ^"nnax  ■  (29) 


From  (28)-(29)  and  p0U 2  =  dax  jds ,  we  obtain 


dT nr 


de 


=  P0V2  +  -4(e)  Tmax  ; 


^  =  (p(p2  -  iPcV2)  -  %A(s)t„ 


A{e)  = 


1  _2  dV 

—  e  V  de 


(30) 

(31) 

(32) 


Since  approximations  for  dV/ds  are  given  by 
(25) ,  for  isotropic  solids  a  single  oblique  plate  im¬ 
pact  test  with  in-material  particle  velocity  gauges 
can  yield  estimates  for  aY ,  rmax,  and  p  in  the  uni¬ 
axially  strained  state  ahead  of  the  shear  wave,  as 
well  as  for  their  derivatives  w.r.t.  uniaxial  strain. 


_  2 e  -  e2  +  72  rxy 
*  v  1  -e  7 

_  2g  ~  g2  Zxy  ,  TJxy 
1  —  c  7  1  —  S 


(26) 


If  the  oblique  plate  impact  generates  a  shear 
shock,  then  evaluating  (26)  in  the  final  state  and 
using  the  jump  conditions  (5)i-(7)i  gives 


(crx  ~  °V)F 


2gF  —  gF2  4-  7f' 

l-sF 

2gr  ~  gr2 

1  —  S'f 


2  ,  Po^r 

1-e. 


(27) 


In  the  weak  shock  limit  7 FSzvf  ->  0,  eF  ->  sl ,  and 
¥  -»  V1 ,  yielding  a  universal  relation  for  the  dif- 
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Abstract.  We  derive  some  approximate  relations  between  shear  wave  speeds  in  inelastic  solids  subject  to 
oblique  plate  impact.  Attention  is  restricted  to  materials  for  which  the  fast  wave  is  purely  longitudinal.  This 
includes  isotropic  solids  as  well  as  appropriately  aligned  orthotropic  solids.  For  the  case  where  the  slower 
wave  is  a  shear  shock,  we  obtain  approximate  relations  between  this  shock  speed  and  the  shear  sound  speeds 
(i.e.,  acceleration  wave  speeds)  immediately  ahead  of  and  behind  the  shock.  These  relations  are  universally 
valid  for  isotropic  or  orthotropic  elastic  solids  as  well  as  inelastic  solids  with  instantaneous  elastic  response. 
Keywords:  Transverse  Waves,  Shear  Shocks,  Acceleration  Waves,  Oblique  Plate  Impact 
PACS:  4640. Cd,  62.50.+p,  81.70.Bt 


INTRODUCTION 

The  waves  studied  in  this  paper  are  of  the  type  gen¬ 
erated  in  an  oblique  plate  impact  test.  A  flyer  plate 
impacts  a  parallel  target  plate,  with  both  plates  in¬ 
clined  to  the  direction  of  motion  of  the  flyer.  The 
particle  velocity  imparted  to  the  target  plate  has  com¬ 
ponents  normal  and  parallel  to  the  impact  surface. 
We  choose  a  Cartesian  coordinate  system  with  the  X\ 
and  X2  axes  in  these  two  directions,  respectively,  and 
restrict  attention  to  points  in  the  target  interior  and 
to  times  prior  to  the  arrival  of  reflected  waves  from 
free  surfaces  or  material  interfaces.  If  the  material 
is  homogeneous  the  impact  generates  plane  waves 
propagating  in  the  X\  direction:  a  quasi-longitudinal 
wave  propagates  into  undeformed  material  at  rest, 
followed  by  a  slower  moving  quasi-trans verse  (i.e., 
shear)  wave.  Here  quasi  refers  to  the  fact  that  the  lon¬ 
gitudinal  and  transverse  waves  are  generally  accom¬ 
panied  by  relatively  smaller  changes  in  transverse 
and  longitudinal  quantities,  respectively.  If  the  mate¬ 
rial  is  isotropic,  however,  the  fast  wave  will  be  purely 
longitudinal.  In  this  case  the  shear  wave  propagates 
into  uniaxially  strained  material,  a  condition  essen¬ 


tial  for  our  main  results.  But  even  for  isotropic  solids 
this  shear  wave  will  still  be  quasi-transverse  due  to 
nonlinearities  in  the  material  response;  cf.  Bland  [1], 
Davison  [2],  Abou-Sayed  &  Clifton  [3],  and  Gupta 
[4]  for  some  of  the  earlier  theoretical  and  experimen¬ 
tal  work  on  this  problem. 

Depending  on  the  materials  and  the  strength  of  the 
impact,  the  waves  generated  could  be  shock  or  ramp 
(i.e.,  continuous)  waves.  We  consider  only  the  case 
where  the  shear  wave  is  a  shock;  whether  the  uniaxi¬ 
ally  strained  state  ahead  of  this  shear  shock  is  gener¬ 
ated  by  a  shock  or  ramp  wave  is  inconsequential. 

Bland  (p.  764  of  [1])  showed  that  to  within  an  er¬ 
ror  of  fourth  order  in  the  jump  in  shear  strain  across 
a  shear  shock,  the  shock  speed  exceeds  the  sound 
speed  ahead  of  it  by  one  third  the  excess  of  the  sound 
speed  behind  it  over  the  sound  speed  ahead  of  it.  Of 
course,  the  sound  speeds  refer  to  shear  waves  as  well. 
Bland  derived  this  result  for  homogeneous  nonlin¬ 
ear  isotropic  hyperelastic  solids,  assuming  the  shear 
shock  is  “propagating  into  an  unstrained  medium 
at  rest.”  But  even  purely  transverse  loadings  on  the 
boundary  would  generate  a  longitudinal  precursor,  so 
his  restriction  is  not  experimentally  achievable. 
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In  [5]  I  noted  that  Bland’s  approximate  relation 
holds  with  no  restriction  on  the  uniaxial  strain  ahead 
of  the  shear  shock,  and  that  “isotropic”  could  be 
weakened  to  “orthotropic  relative  to  the  coordinate 
axes.”  The  proof  (which  was  omitted)  utilized  results 
in  an  earlier  paper  [6],  making  substantial  use  of  the 
hyperelasticity  assumption  (i.e.,  the  internal  energy 
is  a  potential  for  the  stress)  and  also  of  the  method  of 
characteristics  for  calculating  the  sound  speeds. 

Here  I  show  that  Bland’s  approximate  relation  also 
holds  for  a  large  class  of  inelastic  solids  with  prop¬ 
erties  that  may  vary  smoothly  in  the  X\  direction. 
The  proof  makes  no  use  of  a  potential  relation  for 
the  stress  tensor.1  The  method  of  characteristics  for 
calculating  sound  speeds  must  be  abandoned  since 
at  this  level  of  generality  the  equations  of  motion  are 
not  quasi-linear  hyperbolic  PDEs.  Instead,  we  use  the 
theory  of  singular  surfaces  since  the  sound  speeds 
may  be  identified  with  acceleration  wave  speeds.2 


INSTANTANEOUS  ELASTIC  RESPONSE 

The  results  in  this  section  are  stated  for  an  arbitrary 
motion  x  =  x(X,f)  relative  to  an  arbitrary  Carte¬ 
sian  coordinate  system:  x  =  {x\,X2,X'})  is  the  posi¬ 
tion  at  time  t  of  the  material  point  initially  at  X  = 
(Xi,X2,X3).The  deformation  gradient  F  has  compo¬ 
nents  Fij  =  dxi/dXj.  Let  T  and  S  =  (detF)TF-7  de¬ 
note  the  Cauchy  and  first  Piola-Kirchhoff  stress  ten¬ 
sors.  The  value  of  the  stress  at  time  t  can  depend  not 
only  on  F(/),  as  for  an  elastic  material,  but  also  on 
the  entire  past  history  of  the  deformation  gradient: 

S(0  =  S(F(f);FO.  (1) 

Here  Fr  is  the  history  of  F  up  to  time  t ,  the  (tensor¬ 
valued)  function  defined  on  the  positive  reals  by 


1  For  the  materials  considered  here  such  a  relation  would  follow 
from  the  usual  form  of  the  2nd  law  of  thermodynamics;  cf.[7]  for 
materials  with  fading  memory  and  [8]  for  viscoplastic  materials. 
For  simplicity,  thermodynamic  effects  are  neglected. 

2  Acceleration  waves  are  propagating  singular  surfaces  across 

which  the  stress,  strain  and  particle  velocity  are  continuous  but 
suffer  jump  discontinuities  in  their  spatial  gradients  and  time 

derivatives;  they  can  be  the  fronts  of  ramp  waves.  Similarly,  we 
use  the  term  shock  wave  only  for  propagating  singular  surfaces 

across  which  the  stress,  strain  and  particle  velocity  suffer  jump 

discontinuities.  The  jump  in  a  quantity  'F  across  a  singular  surface 
is  denoted  by  [T]  =  —  'F+,  where  T1-  and  T+  denote  limits 

of  from  behind  and  ahead  of  the  wave  front,  respectively. 


F*  (s)  =  F(t  —  s)  for  s  >  0.  Thus  S  is  a  function  of  F (t) 
and  a  functional  of  Fr.3  For  any  orthogonal  tensor  Q 
and  any  tensor  H  in  the  material  symmetry  group,4 

S(QF(f)H;QF  H)  =  QS(F(f);F*)H.  (2) 

We  assume  that  S  is  continuous  at  any  instant  t  at 
which  F  is  continuous.  Across  a  shock  at  time  t ,  we 
assume  that 

S"  =  S(F";F*) ,  S+  =  S(F+;F') ,  (3) 

which  requires  a  continuous  dependence  on  the  past 
history  across  the  shock  front.  At  any  instant  t  when 
F  and  its  material  time  derivative  F  are  continuous, 
we  assume  that  S  is  continuous  as  well  and  given  by 

S  (t)  =  aiS(F(f);F0F(f)  +  G(F(O;F0 .  (4) 

The  4th-order  tensor  d\  S  is  the  derivative  of  S  with 
respect  to  its  first  argument  F(t),  holding  Ff  fixed. 
The  functions  t  i— >  d\S(F{t);Ft)  and  t  ^  G(F{t);Ft) 
are  assumed  to  be  continuous  at  any  instant  at  which 
F  is  continuous,  so  that  across  an  acceleration  wave 

[SJ  =  <?iS(F(f);Ff)  [F] ,  (5) 

analogous  to  the  case  for  an  elastic  material  (where 
the  dependence  on  Ff  would  be  absent).  Materials 
with  the  properties  described  above  will  be  said  to 
have  instantaneous  elastic  response.5 * 

Let  a  =  x  denote  the  particle  acceleration.  Across 
an  acceleration  wave,  balance  of  momentum  and 
kinematic  compatibility  yield  the  relations 

[S]n=  —  p0c|[a],  c[F]4V[a]®n,  (6) 

where  n  is  the  unit  normal  to  the  referential  wave 
front,  c  the  referential  (or  Lagrangean)  wave  speed, 
and  p0  the  initial  density.  On  applying  (5)  to  n  and 
using  (6),  we  obtain  the  propagation  condition 

Poc2H  =  A(n)M-  (7) 

For  a  given  direction  of  propagation  n,  the  2nd-order 


3  The  results  in  this  section  are  unaltered  if  we  also  allow  the 
material  properties  to  depend  smoothly  on  the  position  X. 

4  Indeed,  T  is  unchanged  and  S  changes  to  SH  when  F  is  re¬ 
placed  by  FH  for  all  times.  And  frame-indifference  implies  that 
T  changes  to  QTQ7  and  S  to  QS  when  F  is  replaced  by  QF  for 
all  times.  On  combining  these  results,  we  obtain  (2) . 

5  These  properties  are  possessed  by  materials  with  fading  memory 
[7]  and  also  by  many  internal  state  variable  models  and  models  for 
viscoplastic  response  [8]. 
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tensor  A(n),  called  the  acoustic  tensor ,  is  defined  by 

A(n)b=  («?iS(F(f);F0 [bigin])  n  (8) 

for  any  vector  b.  Thus  (7)  states  that  [[a]]  must  be  an 
eigenvector  of  the  acoustic  tensor  and  that  p0c 2  is  its 
eigenvalue.  The  ij  component  of  A(n)  is  given  by 
Hi  i=\(dSik/ dFji)nkiti,  holding  fixed.  When  n  is 
the  unit  vector  ei  along  the  X\  axis,  this  reduces  to 


A(ei)  = 


dS  11 

dS  11 

dS  11 

dFn 

dFn 

dFn 

dS2i 

dSn 

dS2\ 

dFn 

dF2\ 

dFn 

dSn 

dSn 

dSn 

dFn 

dFn 

dFn  . 

(9) 


PLANE  WAVES 


Assume  the  displacements  depend  on  X\  and  t  only: 


xi-Xi=xi(X1,t),  i=  1,2,3.  (10) 


Then  the  same  holds  for  F,  which  reduces  to 


F  (Xht) 


Fn  0  0 
F2 1  1  0 

F3 1  0  1 


(11) 


For  the  constitutive  relation  (1) ,  S  also  depends  only 
on  X\  and  t.  For  fixed  X\  we  may  write  this  as 

Sij (0  =  SijiFi  1  (?) ,  Fn  (?) , F31  (?) ;  F(x ,  F{x ,  F{x ) . 

(12) 

Momentum  balance  reduces  to 


we  see  that  Sij  =  —Sij  when  F3\  =0.  Thus  S  reduces 


to 


S  (Xut) 


s  11  sl2  0 

^21  S22  0 

0  0  S33 


(15) 


when  F31  =0.  In  this  case  we  also  have 


dS3\ 


dFn 


f31=o 


dS3\ 

dF2i 


=  0. 

^31  =  0 


(16) 


Now  assume  that  V3  =  X3,  which  is  consistent  with 
the  boundary  conditions  discussed  in  the  Introduc¬ 
tion.  Then  F31  =  0  and  a3  =  0,  so  the  3rd  component 
of  the  momentum  balance  equation  (13)  is  satisfied 
since  S3\  =0.  By  (14)  we  see  that  S n  and  S2\  are 
even  functions  of  F3\  when  Fix  =  0.  Thus  we  also 
have 


dS  11 


dF3 1 


^31  =  0 


dS2\ 
dF3 1 


=  0. 

^31  =  0 


(17) 


Define  the  longitudinal  stress  and  strain,  a  and  £,  and 
the  shear  stress  and  strain,  t  and  y,  by 


a  =  —S 11  =  -Fn 
e  =  1 


dxi  _ ,  F 

wr  111 


T 

r 


-S21  = 
dx2 
~dx~l 


-721, 

(18) 

—7^1- 


These  variables  should  be  nonnegative  for  the  impact 
problem  considered  here,  and  we  may  now  write 


{<?1}  =  {?}(£W’rW;£’r?f  (19) 

By  (16) -(18),  the  acoustic  tensor  (9)  for  waves 
propagating  in  the  X\  direction  reduces  to 


/=  1,2,3,  (13) 

which  involves  only  the  stress  components  S\\,  S2\ 
and  S3\ .  This  continues  to  hold  if  the  material  prop¬ 
erties  vary  smoothly  with  X\ .  The  motion  is  a  plane 
wave  propagating  in  the  X\  direction. 

Now  assume  the  symmetry  group  contains  a  180° 
rotation  about  the  X3  axis;  its  matrix  is  given  by 
diag(— 1,  —  1, 1).  On  choosing  this  for  Q  and  H  in 
(2) ,  we  obtain  the  following  restrictions  on  (12) : 


A(ei)  = 


da 

da 

~de 

dy 

dT 

dT 

de 

dy 

0 

0 

(20) 


0 

0 

<?S3i 


dF31 


£31  =  0. 


Since,  by  (7),  p0c2  is  an  eigenvalue  of  A(ei),  two  of 
the  Lagrangean  sound  (acceleration  wave)  speeds  are 
obtained  by  equating  2 p0c2  with 


Sij(FiuF2u-Fn-,F{l,F^,-F^)  (14) 

=  (En,  7*2!,  E31;  7*n ,  7*21 , 7*3^ )  > 

where  the  +  case  holds  for  ij  =  11, 22, 33, 12, 21,  and 
the  —  case  for  ij0  13,23,32,31.  For  the  latter  case 


da 

~de 


dr  (3°  dT\2  ddo  dT 

dY±y\fc~dy)  +4lfyd£' 


(21) 


The  faster  (quasi-longitudinal)  sound  speed  is  de¬ 
noted  by  U,  and  the  slower  (quasi-trans verse)  sound 
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speed  by  V.6 7 

Now  assume  the  symmetry  group  also  contains  a 
180°  rotation  about  the  X\  axis,  so  that  the  material 
is  orthotropic.  On  taking  Q  =  H  =  diag(l,  —  1,  —1) 
in  (2),  we  obtain  restrictions  on  Sij  in  (12)  which, 
in  view  of  F^\  =0  and  (18)— (19) ,  imply  that  for  any 
strain  values  £*  and  y*  and  strain  histories  and  yl, 

9(£*,-r*;e|,-7')  =  S(£’,r-;£',r'); 

It  follows  that  if  y  =  0  up  to  and  including  time  t , 
then  so  are  x,  dx/de  and  da /dy? 


SHEAR  SHOCKS 

Consider  a  shear  shock  propagating  into  uniaxially 
strained  material  behind  the  longitudinal  wave  in 
an  oblique  plate  impact  test.  Let  ¥  denote  the  La- 
grangean  shear  shock  speed.  The  jump  conditions  are 

W=P0V2W,  M=PoV2M;  (23) 

hence 

Wll£]  =  M[Irl-  (24) 

Since  y  and  x  are  zero  ahead  of  the  shock,  we  have 
y+  =  T+  =  0.  Thus  [[y]]  =  y~ ,  and  by  (22)  the  longi¬ 
tudinal  and  shear  stress  behind  the  shock  are 

cr“=a(e“,[[7];£f,0), 

M  =  t-=?(£-,[[y]]  ;£',()), 

these  being  even  and  odd  functions,  respectively,  of 
[y]  for  fixed  e~  and  .  But  e~  and  [[y]]  are  coupled 
through  (24) -(25),  which  can  be  used  to  show  that 
[e]  is  of  order  [[y]]2.  Then  (23)  i  and  the  sound  speed 
relations  imply  ¥  ~  V+  ~  V-  to  within  an  error  of 
order  [[y]]2,  where  V_  and  V+  are  the  limiting  values 
of  the  shear  sound  speeds  from  behind  and  ahead  of 


6  Note  that  these  sound  speeds  are  well-defined  at  each  point  and 
time,  regardless  of  whether  or  not  the  corresponding  acceleration 
wave  is  actually  generated  in  the  impact  test.  The  third  speed, 
given  by  p0c 2  =  (dS3i/dF2i)\F3l=o,  corresponds  to  a  shear  wave 
polarized  along  the  X3  axis  and  is  not  relevant  here. 

7  Since  £21  =  —  T,  the  fact  that  T  =  0  for  a  pure  uniaxial  strain 
history  suffices  for  the  fast  wave  in  an  oblique  plate  impact  test 
to  be  purely  longitudinal.  The  zero  derivatives  imply  that  the 
relations  for  the  longitudinal  and  shear  sound  speeds  in  uniaxially 
strained  material  reduce  to  p0U 2  =  do /de  and  p0V2  =  dx/dy, 
evaluated  at  y  =  0;  since  (20)  is  diagonal,  these  are  pure  modes. 


the  shear  shock.  It  follows  that  ¥  «  aV+  +  /3  V_  for 
any  a  and  /3  satisfying  a  +  /3  =  1 ,  also  to  within  an 
error  of  order  [[y]]2. 

It  turns  out  there  is  an  optimal  choice  for  a  and  /3 
above.  From  (21)  and  the  results  in  this  section,  it  is 
not  hard  to  show  that  to  an  within  error  of  order  [[y]4, 


2  ^  dx  da  dx  I  f  da  dx\ 

~  dy  dy  de  /  \  de  dyj  ’ 


(26) 


where  the  right  side  is  evaluated  behind  the  shock, 
at  (e-,  [[y]];e*,0).  Some  lengthy  calculations  reveal 
that  this  equals  p0(3¥2  —  2¥+)  to  within  an  error  of 
order  [[y]]4.  Hence,  we  obtain  the  equivalent  relations 
below,  all  valid  to  within  an  error  of  order  [y]4: 

¥2  «  |V2  +  2y2  ¥2-V2  «  \{Vl-Vl) 

(27) 

V2  ^3¥2-2V2  V2-¥2  ^2(¥2-V2). 


These  yield  analogous  approximations  for  the  speeds 
themselves  to  within  an  error  of  order  [y]4: 

¥ «  |V+  +  \V-  ¥-y+^  2(y__y+) 

(28) 

y_  ^3¥-2y+  y_-¥^2(¥-y+) 


These  results  depend  in  an  essential  way  on  the 
condition  of  uniaxial  strain  ahead  of  the  shear  shock. 
Although  we  have  neglected  thermodynamic  effects, 
it  can  be  shown  that  their  inclusion  does  not  alter  the 
above  results  if  we  make  the  typical  assumption  that 
the  heat  flux  does  not  jump  across  a  shock. 
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D  CASEM 
J  CLAYTON 
D  DANDEKAR 
M  GREENFIELD 
Y  HUANG 
BLOVE 

M  RAFTENBERG 
E  RAPACKI 
M  SCHEIDLER 
T  WEERASOORIYA 
AMSRD  ARL  WM  TE 
B  RINGERS 
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